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Four-pion production* 

G. Ecker and R. Unterdorfer a 

a Institut fur Theoretische Physik, Universitat Wien 
Boltzmanngasse 5, A-1090 Vienna, Austria 

Starting from the low-energy structure derived from QCD, we extend the amplitudes for four-pion production 
in e + e~ annihilation and r decays up to invariant four-pion masses of 1 GeV. Cross sections and branching ratios 
BR(p° — > 47r) are compared with available data. 



1. INTRODUCTION 

The production of four pions in e + e~ annihila- 
tion and in r decays is interesting in its own right 
but it also represents a non-negligible component 
of hadronic vacuum polarization. In fact, almost 
5 % of the lowest-order hadronic contribution to 
the anomalous magnetic moment of the muon 
is due to four pions 0. At the level of accuracy 
required for a comparison of the standard model 
prediction for with the recent BNL measure- 
ment ||, better knowledge of this contribution 
would be welcome (see the discussion in Ref . ) . 
For the determination of the fine structure con- 
stant at s = Af|, the relative importance of the 
four-pion contribution is even bigger but in this 
case the available precision is sufficient for the 
time being. 

In this talk, I report on the work of Ref. || 
where we have constructed the relevant ampli- 
tudes up to invariant four-pion masses of about 1 
GeV. This is not enough for the purpose of cal- 
culating a M , but it is a first step in this direction. 
Even if it is impossible to calculate the amplitudes 
directly from QCD our aim was to construct am- 
plitudes that are at least consistent with QCD. 
This program turned out to have some surprises 
in store. 

The procedure of Ref. || starts from a calcu- 
lation to next-to-leading order in the low-energy 
expansion of QCD, employing the methods of chi- 
ral perturbation theory (CHPT) 0]. The corre- 
sponding low-energy amplitudes cannot be used 



directly in the physical region because the four- 
pion threshold of 560 MeV is already close to the 
resonance region. But the low-energy amplitudes 
contain nontrivial information how to continue 
to higher energies. At next-to-leading order, the 
traces of p and scalar meson exchange appear in 
the amplitudes. Supplemented by to, a\ and dou- 
ble p exchange, the resulting e + e~ cross sections 
describe the available experimental data very well 
up to cms energies of about 1 GeV. 

2. SYMMETRIES AND LOW-ENERGY 
LIMIT 

There are altogether four different channels ac- 
cessible in e + e~ annihilation and r decays into 
four pions. In the isospin limit that we assume 
throughout, the amplitudes of either e + e~ — > 
27r°7r+7r- [§ or t~ -> v T 27r-7r+7r° [§ arc suf- 
ficient to determine all four amplitudes. One im- 
portant advantage of the chiral approach is that 
not only chiral symmetry but also charge conjuga- 
tion invariance, Bose symmetry and electromag- 
netic gauge invariance are manifest at each stage 
of the calculation. 

At leading order in the low-energy expansion, 
the amplitude is completely determined by "vir- 
tual" bremsstrahlung. From the diagrams in 
Fig. ^, the matrix element of the electromag- 
netic current governing the amplitude for e + e~ — > 
27r°7r + 7r~ is found to be 
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Figure 1. Tree diagrams for 7* — > in. Solid lines 
denote pions and the wavy line stands for the vir- 
tual photon. 



where s = (pi + p 2 ) 2 , q — Pi + P2 + P3 + Pa and 
= 92.4 MeV is the pion decay constant. This 
matrix element has a very suggestive structure: 
(s — M%)/F% is the (lowest-order) scattering am- 
plitude for 7T°7r — * 7T+7T - and the second term re- 
duces to the usual bremsstrahlung factor for real 
photons (q 2 — > 0). 

Although the amplitude (0) is by itself not of 
direct phenomenological relevance we can use the 
isospin relations || and calculate also the leading- 
order t decay amplitudes. In the chiral limit 
(M„ = 0), those amplitudes should coincide with 
those of Ref. |J. Unfortunately, there are two 
(identical) misprints in Ref. || that have more- 
over propagated into some of the subsequent lit- 
erature (0 and references therein). Although 
the structure is correct, the normalization is not: 
the correct amplitude (0) and the correspond- 
ing t decay amplitudes are smaller by a factor 
\/2/(3v3) (or 1/13.5 in rate). This normaliza- 
tion error also affects the so-called CLEO current 
f§ in the Monte Carlo package TAUOLA @. 

The "current algebra" amplitude (Q) is impor- 
tant for checking the low-energy limit of QCD but 
it is not a realistic approximation in the physical 
region. In order to see the traces of meson reso- 
nance exchange, we have to go at least to next- 
to-leading order. 

3. RESONANCE EXCHANGE 

At next-to-leading order in the chiral expan- 
sion, the amplitude consists of two parts: a loop 
amplitude JlC(] and a tree amplitude containing 
the (renormalized) coupling constants of the chi- 



ral Lagrangian of 0(p 4 ) Q. With the standard 
values of those constants, one arrives at cross sec- 
tions that are still unrealistic. As indicated by the 
dotted curves in Figs. 2,3, the theoretical cross 
sections are significantly smaller than the mea- 
sured ones. 

The seeds of meson resonance exchange ap- 
pear first in the coupling constants of 0(p A ). In 
fact, those constants are known to be saturated 
by meson resonance exchange to a large extent 
fi"l] , [r2|| . This saturation makes the matching be- 
tween the strictly chiral amplitude to 0(p 4 ) and 
a more realistic meson resonance exchange am- 
plitude almost trivial. Using the standard chiral 
resonance Lagrangian pd[ , the resonance ampli- 
tudes are guaranteed to exhibit the correct low- 
energy behaviour to 0(p 4 ). 

In four-pion production, only the p and the 
(isoscalar) scalar mesons contribute at 0(p 4 ). As 
could have been expected, p exchange dominates 
by far. The overall contribution from scalar ex- 
change turns out to be very small so that the 
controversial structure of the scalar sector is not 
relevant in practice. 

The modified amplitudes with p and scalar ex- 
change are definitely more realistic than the chi- 
ral low-energy amplitudes. However, except in 
the vicinity of the p pole, the resulting cross sec- 
tions are still too small (not shown in Figs. 2,3). 
The obvious lesson is that important ingredients 
of the amplitudes are still missing that only show 
up at orders p e or higher in the chiral expansion. 

The most important missing degrees of free- 
dom are easily found: both quantum number con- 
siderations and experimental information 13 . 14] 



indicate that oj and a\ exchange must be incor- 
porated. Whereas the lowest-order coupling of 
the u> to pions is unique, there is some ambiguity 
in the aipn couplings (to be resolved eventually 
by studies of three-pion production |15|). With 
a simplifying assumption for those couplings ||, 
including double p exchange that also comes in 
at 0(p 6 ) and performing a resummation of some 
terms making up the p-dominated pion form fac- 
tor, the amplitudes assume their final form. Ex- 
cept for the ambiguity in the aipn couplings, all 
resonance couplings can be determined from the 
respective decay widths. 
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4. COMPARISON WITH DATA 

The two main assets of our amplitudes are: 

• They contain the relevant degrees of free- 
dom for describing four-pion production up 
to energies of about 1 GeV. 

• They exhibit the correct low-energy be- 
haviour to 0(p 4 ) by construction. 

For energies below 1 GeV, annihilation data 
are available for the channel 27r + 27r~ mainly. In 
Fig. 2, the theoretical cross sections are compared 
with the most recent (and most precise) data from 
the CMD-2 Collaboration jw| (see Ref. § for the 
full data set). The cross section for our model is 
shown as the full curve. The dashed curve corre- 
sponds to omitting the loop amplitude of 0(p 4 ) 
|l0[ (except for the contribution to the width of 
the p meson) . The obvious conclusion is that the 
amplitude is completely dominated by resonance 
exchange. Although a possible enhancement of 
the cross section in the region between 800 and 
900 MeV could not be explained with our am- 
plitude the gross features of the data can be re- 
produced over a range of two orders of magni- 
tude with almost no free parameters. In retro- 
spect, the simplifying assumption for the aipir 
couplings (that actually has a theoretical basis 
15 1) is justified by comparison with experiment: 
other choices for the couplings could not repro- 
duce the data. Note that w exchange does not 
contribute in this channel: the cross section near 
1 GeV is completely dominated by a\ exchange. 

The experimental situation is less satisfactory 
for the other annihilation channel 27r 7r + 7r~. The 
most precise experiment |L3j has measured the 
cross section at only two energies below 1.05 GeV. 
The comparison between theory and experiment 
is shown in Fig. 3. The loop contribution is 
even less relevant in this case. The theoretical 
cross section increases by three orders of mag- 
nitude from the p resonance to match the two 
data points. For this channel, the cross section 
near 1 GeV is dominated by u exchange. The 
theoretical prediction is therefore less sensitive to 
assumptions about the aipir couplings. 

Our amplitudes and the corresponding cross 
sections cannot be extended to the phenomeno- 



logically most interesting region above 1 GeV 
without further input. The reason is that the am- 
plitudes do not satisfy the high-energy constraints 
of QCD. In fact, the theoretical cross sections ex- 
ceed the data soon above 1 GeV of cms energy. 
Additional higher-mass states must be included 
to access the region up to 2 GeV and to ensure 
a proper high-energy behaviour. Rcsummations 
similar to the pion form factor may also be nec- 
essary. 




0.8 0.85 0.9 
energy (GeV) 



Figure 2. Comparison of theoretical predictions 
for e + e~ — > 2tt + 2it~ with data [jl4|: complete 
cross section (full curve), without loop contribu- 
tions (dashed curve), cross section for the ampli- 
tude of 0(p 4 ) (dotted curve). 



From our amplitudes we can also extract the 
branching ratios for the four-pion decay modes of 
the p°. For q 2 = M 2 several contributions to the 
amplitudes are of comparable size, with partly 
destructive interference. In this way, uncertain- 
ties in the resonance couplings arc enhanced. We 
therefore quote predictions for the branching ra- 
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Figure 3. Same as in Fig. 2 for the process 
e+e- -> 27r°7T+7r- p|. 



tios with a 40 % uncertainty: 



BR(p° -> 2tt+2tt-) 
BR(p° -> 27r°7r + 7r") 



(6.7 ±2.7) x 10~ 6 (2) 
(5.0 ±2.0) x 10~ 6 . (3) 



For comparison, the Particle Data Group [|16| lists 

BR(p° -> 27r+27r-) = (1.8 ± 0.9) x 10~ 5 (4) 
BR{p° -> 27r°7r + 7T") < 4 x 10~ 5 . (5) 

5. CONCLUSIONS 

The following features of our model for four- 
pion production are worth repeating: 

• The amplitudes exhibit the correct low- 
energy structure to 0(p 4 ) in the chiral ex- 
pansion. 

• All symmetries of the transitions are man- 
ifest in the QFT framework of CHPT: 
(broken) chiral symmetry, gauge invariance, 
Bose symmetry and charge conjugation. 



• In addition to p (and the less important 
scalar) exchange, to and a\ exchange are 
crucial for understanding the experimental 
results already at energies below 1 GeV. 

• Good agreement with available data is ob- 
tained, covering several orders of magnitude 
in cross sections. 

• To extend the amplitudes to energies above 
1 GeV, the correct high-energy behaviour 
still needs to be implemented. For the same 
reason, comparison with r decay data is 
postponed. 
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